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Abstract 

Abstract: Under rather general assumptions, and in a relatively 
simple and straightforward manner, it is shown that the characteris- 
tics of signals which travel through homogeneous, as well as inhomo- 
geneous, passive media have the principal features usually associated 
with the phenomena of precursors, as generally follows from more elab- 
orate studies. The simplicity of the present arguments permit analytic 
studies to be made for a greater variety of media than is normally the 
case. 



1. Introduction 

The propagation of electromagnetic signals through passive media, i.e., media 
exhibiting both dispersive and absorptive properties, of many different kinds 
is a problem of interest in many applications. For simplicity of analysis, we 
model such problems in relatively economical but fairly general terms. Let 
fo(t) denote the temporal dependence of the initial signal at some initial 
depth level (say, z = 0) inside, or at the surface, of some medium. We 
presume the medium to be homogeneous at this stage and that its response 
characteristics are stationary. We denote by f z (t) the temporal dependence of 
the signal after traveling a distance z > through the passive medium. The 
relation between the input and output signals is that of a standard nitration 
which may be presented in several different forms. One such form is given 
by convolution with the medium impulse response function m z (t), and takes 
the form 

/*(*)= / m z (t-s)f (s)ds. (1) 
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Note that all integrals with unspecified limits are assumed to extend over 
the whole real line. As a convolution, these functions have a multiplicative 
relation in frequency space given by 

F g (u) = M z (uj)F (uj) , (2) 

where 

FM = f e^f z (t)dt, (3) 

M x (u) = J e™m g (t)dt, (4) 

denote the corresponding Fourier transforms. As a causal phenomena (no 
signal out before a signal in), it follows that m z (t) = 0, for t < 0, a property 
that means that M z (u) has no poles in the upper half complex u plane, i.e., 
no poles for oo = uo r + iuoi with uoi > 0. These few properties are standard for 
causal signal transmission [1]. 

Let us now take up an elementary property of M z {uj) regarding its de- 
pendence on the depth parameter z. Clearly, 

F Zl+Z2 {u) = M Zl {u)F Z2 {uo) 

= M Zi {uj)M Z2 {uj)F q {uj) 

= M Zi+Z2 (u)F (lu) . (5) 

Since F (cu) is arbitrary in this expression, it follows that 

M Zi+Z2 (lo) = M Zi (uj)M Z2 (lo) (6) 

holds for all nonnegative z\ and z 2 . In particular, it follows that 

MqH = MoH 2 (7) 

implying that either Mq(uj) = 1 or Mq(uj) = 0. If there are frequency bands 
that F (uj) will never possess, then it is possible to choose Mq(u>) = in that 
band, but for simplicity we may as well assume that M (u) = 1 for all ui 
since the final answer will be unchanged. Moreover, on physical grounds it 
is natural to assume, for each u>, that M z {uj) is a continuous function of z. 
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In that standard argument asserts that the only function that fulfills 

(6) necessarily has the functional form 

M z {uj) = e~ zA ^ . (8) 

Let us recall that argument. 

To establish the stated form for M z (u), we first observe that 

M z (u) = [M g/n (u)] n (9) 

for an arbitrary positive integer n. Since M z {uj) is assumed to be continuous 
and Mq{u) = 1, it follows that M z i n {uj) ^ for large enough n, and so 
M z (uj) 7^ for any z. Next, let R z (u) = hi(M z (u)), which is well defined for 
all z and uj. Furthermore, (6) leads to R Zl+Z2 (u) = R zi {uj) + R Z2 (cu). This 
relation implies that R nz {uS) = nR z (u>), as well as R z / m {uj) = (l/m)R z (u), 
for arbitrary positive integers n and m. Consequently, 

RnzimW) = (n/m)R z (uj) . (10) 

Now let z — 1 so that R n / m {oj) = (n/m)Ri(u); then take a limit such that 
(n/m) — > z, an arbitrary positive number. By the assumption of continuity, 
this limit converges and we learn that 

R z (ou) = zRi(w) = -zA(w) . (11) 

This concludes the proof of (8). 

As the transfer function of a passive medium it is necessary that \M z (uo) \ < 
1, in which case we also find that JieA(uj) > 0. To simplify notation, we let 

A(u) = B(lo) + iC(uj) , (12) 

where B(u) = ReA(u) and C{uo) = \m.A(uo). Thus, we require that 

B(u) > . (13) 

By causality, B{uj) and C{uj) are related by Hilbert transforms, while C(u) is 
undetermined by this procedure up to the contribution of a so-called Blaschke 
factor [2]. We also note that, for real functions f z (t) for all z and any initial 
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signal /o(t), it necessarily follows that the impulse response function m z (t) 
is a real function. Consequently, M 2 (— u) = M z (cu)*, which implies that 

B(-u) = B{u) , (14) 
C(-u) = -C{uj) . (15) 

Since B(u>) > 0, it follows that the function B{uj) has special properties at 
and near any frequency, say u/, for which B{uj') = 0. Assuming on physical 
grounds that B{uj) is continuous in u, and even has a few derivatives [e.g., 
B(uo) E C 2 ], then the behavior of B in the vicinity of oo' is given by 

B(u) = \a-\u-u') 2 , (w-a/) 2 <l. (16) 

where a > has the units of acceleration (e.g., m/s 2 ). 

Let us focus on the behavior of A(u) in the vicinity of to — 0. In that 
case, and to sufficient accuracy, we have 

B(uj) =r 1 + \ar 1 u 2 , (17) 
C(w) = , (18) 

where £, v, and a denote constants with the dimensions of length, velocity, 
and acceleration, respectively. The behavior at uo = forces I" 1 > 0. More- 
over, when B(0) = l~ x = 0, it follows that a -1 > 0; in that case, we shall 
adopt the generic form and assume that a > 0. 

It is natural to ask the question: Under what circumstances does it follow 
that B(0) = 0? A common model of a passive medium is based on the 
Lorentz model [3]. The physical principle behind this model is the response 
of a damped harmonic oscillator to a driving force. The physical model in 
this case is given by 

mq(t) + bq(t) + kq(t) = h(t) , (19) 

where m > 0, b > 0, and k > denote the inertial, damping, and restora- 
tive parameters, q(t) denotes the oscillator amplitude, and h(t) denotes 
the applied force. Assuming that any transients have decayed away, let 
us focus on the steady state solution. If h(t) = U cos(cut), then q(t) = 
V cos(ujt) + W sin(ujt), where 

v = ( k ~ muj2 ) u m 

{k - muo 2 ) 2 + {bu) 2 ' 1 J 

W = ^ . (21) 

(k - moo 2 ) 2 + (bco) 2 y } 
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If the driving frequency uj = for this model it follows that V = U/k and 
W = 0. Since this result is entirely independent of the dissipation coefficient 
b, it holds even when b = 0, namely, it holds in the case there is no loss 
whatsoever. Consequently, for such a model it follows that at uj = there 
can be no loss and so B(0) = 0. Of course, such an argument does not apply 
to a metal, but it does hold for a broad spectrum of materials. 

Although we have argued that B(0) = for a specialized model, it is 
quite plausible that such a result holds for a wide class of models as well. 
One such more general model would be of the form 



which corresponds to a situation where h(t) acts as the driving force on a 
collection of coupled damped harmonic oscillators. Again we may consider 
the steady state response when h(t) = U cos (cut). In particular, if uj = the 
solution is totally independent of any dissipation parameters. Consequently, 
it follows that 5(0) = in this case as well. It is not difficult to imagine other, 
more complicated situations for which the same result holds. Therefore, on 
the basis of a wide class of physically plausible models, we are encouraged 
to consider the consequences of the assumption that B(u>) = |a _1 u; 2 , a > 0, 
when uj is small. 

We now examine the signal behavior as a function of z that follows from 
the foregoing discussion. In particular, we first focus on the impulse response 
function 



where A(oj) = B(u)+iC(u), and B(oj) = 0(uj 4 ) while C(oj) = 0{uj 3 ). To 
simplify matters, let us further assume that B{uj) is bounded away from for 
nonzero values of uj. In that case, for large z, we observe that the integrand 
tends to be strongly suppressed for ui away from zero. To sufficient accuracy 
it follows that we may set 



rriiQi(t) + hqiit) + S J=1 % ^(f) = h(t) , 



1 < % < J , 



(22) 




(23) 




(24) 
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Consequently, for large z values we learn that 



f z (t)= I m z (t-s)f (s)ds 

/ exp { — J fo(s) ds 

J e~ as2/2z f (s + t - v^z) ds . (25) 



2-7TZ 



2nz 

There are two characteristic features of this resultant behavior that are 
worth emphasizing: 

1. According to the square root prefactor, it follows that the amplitude 
decays proportional to z^ 1 ^ 2 . 

2. According to the integral factor, the result has a temporal depen- 
dence determined by the form of the original signal centered at a temporal 
value of t— v and then integrated over that region weighted by a Gaussian 
of increasing width. 



1.1 Approximate causality 

Before passing to several examples of (25) it is important to observe that 
the medium impulse response function (24) is not causal since it does not 
vanish for t < 0. In point of fact, the expressions (17) and (18) do not 
make a strictly causal impulse response function. However, they make an 
approximate causal function provided that az ^> v 2 since in that case the 
impulse response function is extremely small for t < 0. Let us examine this 
situation more closely. 

In order for the medium function A(u) to correspond to a causal function 
it should be represented in the form 

POD 

A(u) = / [1 -e ws ]k(s)ds , (26) 
Jo 

where A(0) = 0, and in order that ReA(u) = B(u) > 0, it is sufficient that 
the weight function k(s) > 0. Observe, by this representation, that A(u>) is 
analytic in the upper half plane, which is a necessary and sufficient condition 
for the impulse response function to be causal. Now, as an example, consider 
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k(s) = K'e Ks , where K' and K are positive parameters to be determined. 
In that case, 

poo 

B(u)=K' [1 -cos(ujs)]e- Ks ds , (27) 
Jo 

poo 

C{uo) = -K' sin(Lus)e- Ks ds . (28) 
Jo 

We are interested in a large value for both K and K' such that, to some 
approximation, k(s) acts rather like a delta function in the integrand focusing 
on values of s near zero. Thus, to a leading approximation, 

roc 

B(uj) = \K'uj 2 / s 2 e~ Ks tfe + • • • , (29) 
Jo 

poo 

C(uj) = -K'uj se- Ks ds + --- . (30) 
Jo 

On comparison with (17) and (18) we see that 

roc 

a' 1 = K 1 / s 2 e~ Ks ds = 2K'/K 3 , (31) 
Jo 

poo 

v- 1 = K' s e~ Ks ds = K'/K 2 , (32) 
Jo 

and so a = \Kv. Provided K ^> 1, this approximation is satisfactory from 
a practical viewpoint; otherwise, some of the additional terms in (29) and 
(30) need to be taken into account. Unfortunately, introducing additional 
terms makes the medium impulse response function no longer analytically 
tractable, generally speaking. 



1.2 Specific examples 

Let us illustrate the resultant signal (25) for two important examples. First 
we consider an analytically simple case, namely a Gaussian initial pulse. 
Thus, let 

foit) = e~ t2 ' 2T2 cos(ua) , (33) 
where we refer to T as the pulse width. In this case, (25) reads 



fz{t) = \IWz 1 e-<°-^-^' 2 *e- s2 ' 2T * cosM ds . (34) 



An elementary integration leads to 
/,(*) ! 



y/l + z/aT 2 

[(a/z)(t 



x exp 



(- 



-M 2 + ^ T 2 



\ fU a (t-V 1 Z)\ 

) cos ( l + ^ )' (35) 



2(1 + aT 2 /z) 



+ z/aT 2 



If we assume that z ^> aT 2 , then in effect (35) reduces to 



fz(t) 



1T 2 



exp 



a(t 



,-i 



zf (u Tf 



2z 



(36) 



The peak amplitude of the resultant pulse f z (t) occurs when z = vt, or, 
stated otherwise when t = v~ x z. At that point in z (or at that moment in 
time t), the amplitude has been reduced from its initial maximum of unity 
by a factor proportional to z~ x l 2 . The resultant pulse shape in the present 
case is still Gaussian, but it is substantially broadened since the pulse width 
T of the input Gaussian has been replaced by \Jzja 3> T. If we assume the 
energy (or power) in a given pulse f z (t) is proportional to J\f z (t)\ 2 dt, then 
it follows, for z 3> aT 2 , that the normalized signal energy is given by 

2,1— : , mn . (37) 



J\f z (t)\ 2 dt 



f\fo(t)\ 2 dt "V z l + e-M 3 

As anticipated, this ratio decreases as z increases. However, thanks to the 
pulse broadening, it decreases only as z^ 1 ^ 2 rather than z -1 , as may have 
been expected on the basis of (25). 

As a second example, we suppose that 

f (t) = rect (t/T) cos(cu t) , (38) 

which leads to the output signal 

r T/2 



hit) = f 2 e^-t^W* cosM ds . 

V Zttz J_ t/2 



(39) 



This integral depends on z in a moderately involved way, and for convenience 
we only evaluate (39) in the case that z ^> aT 2 . Then we find 

r-T/2 



hit) 



2nz 



e -a{t-v^zf/2z CQsMds 

-T/2 



aT 2 

2ttz 



■T/2 

sin(u; T/2) 
KT/2) 



3 -a(t— y" 1 *) 2 / 22 



(40) 
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Again, the result is that of a traveling pulse moving essentially at the speed 
v, the amplitude of which falls off with distance like z -1 / 2 , provided of course 
that sin(cj Q T/2) ^ 0. For the case that sin(cj D T/2) = 0, u T ^ 0, see below. 

1.3 More general inputs 

Apart from an unimportant coefficient, observe that the amplitude factor 
arising from fo(t) in (36), i.e., exp(— u 2 T 2 /2), is basically -F (0) = J fo(s) ds, 
and the amplitude factor in (40), i.e., sin(a> T/2)/(a; T/2) is essentially just 
Fq(0) = f fo(s)ds. We can understand and generalize this result as follows. 
Since 



it follows, for an initial signal which is much narrower in the temporal do- 
main than the medium response function, that we can make a Taylor series 
expansion of the medium response function. In particular, we have 



/*(*) = f[m z (t)f (s) - m' z (t)sf (s) + \m" z {t) s 2 fo{s) + ---}ds 

= m z (t) J f (s)ds - m' z (t) J sf (s)ds + \m" z {t) J s 2 f (s)ds H . 

(42) 



When m z (t) is sufficiently broad in time, then its time derivatives m' z (t), m"(t), 
etc., are typically much smaller. In this case, the indicated sum involves terms 
that are rapidly diminishing relative to the former terms. 

As an example of the application of (42), suppose that uo Q T = 2nir, 
n G {1,2,3, . . .}, then (40) vanishes because J fo(s)ds = 0. Moreover, since 
fo(t) in (38) is an even function, it also follows that Jsfo(s)ds = 0. Thus, 
according to (42), we find to leading order that 



f z (t) = fm z (t- s)fo(s) ds , 



(41) 



fM = \<{t)js 2 h{s)ds, 



(43) 



namely, 



/.(*) 



2(mr) 



-1 




— z e 



-a(t-z/v) 2 /2z 



(44) 



where the amplitude has been evaluated when oo T = 2mr. 
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1.4 Commentary 



Observe, that although the original Gaussian and the rectangular envelopes 
are fundamentally different, the pulses that result after a long propagation 
distance z ^> aT 2 are essentially the same, differing only in their respective 
amplitudes. Those amplitudes effectively measure the spectral content of 
the original pulse at zero frequency where the narrow window of transmis- 
sion exists, and for the Gaussian envelope that content is much smaller than 
is available with the rectangular envelope. This result is entirely consistent 
with the concept that the output selects its preferred signal shape consistent 
with the frequency window allowed by the medium. Ultimately, this means 
that the transmitted signal is Gaussian in character with an ever increasing 
pulse width in accord with the ever narrowing Gaussian pass band as z in- 
creases. This kind of temporal behavior is essentially the kind of behavior 
typically ascribed to the so-called Brillouin precursor [3]. It is paradoxical, 
but true, that Brillouin precursors require pulse envelopes with high fre- 
quency content. Sharp edges in pulse envelopes significantly enhance the 
high frequency content as compared to smooth envelopes, and this situation 
is clearly illustrated in Eqs. (36) and (40). 

However, sharp edges in pulse envelopes are not the only way to raise the 
high frequency content. 

2. Chirp signals 

As an initial wave form consider the linearly chirped signal [4] 



where the signal envelope exit) has a pulse width of T and we assume that 
aT 2 3> I. The Fourier transform of (45) may be approximately given by 



f (t) = e T (t) cos(uj t + \at 2 ) , 



(45) 




(46) 



The frequency content at uo = is then given by 



F (0) = -J-e T (-Lu o /a)e^ 2a + ±- 
2 V %ol 2 




2ixi 



a 



e T (-u /a)e 



-iwl/2a 



(47) 
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As an illustration, consider exit) = exp(— t 2 /2T 2 ). In this case, 

F T (woT) 2 i r f(uJoTf\ /(^ T) 2 \i 
F„(0) = ^- exp [ - yj,] [ cos (L^L) + sm (L^)] . (48) 

As noted above, this factor essentially dictates the amplitude for the output 
Gaussian in (24). 

The main point to observe in this calculation is the significant enhance- 
ment of the overall amplitude of the signal that will emerge. In particular, 
the primary change has been that the factor exp[— (uo T) 2 /2] in (36) has been 
replaced by exp[— {uj T) 2 /2a 2 T A ], where aT 2 ^> 1. This model calculation 
strongly suggests that combining chirp signals with envelopes having sharp or 
nearly sharp edges should lead to signals with relatively significant precursor 
content. 

3. Composite media 

Up to this point we have assumed that the signal propagates through a 
single, homogeneous medium, the properties of which are characterized by the 
complex medium function A(uj). Suppose instead we deal with a composite 
media made up from J layers, each of thickness h — z±, I2 — z-i — z±, . . . , lj — 
zj — zj-\. To each layer we attribute a medium function Aj(u>), 1 < j < J. 
If we assume that Zj < z < Zj+i, for example, then 

M z {u) = exp[-(z - Zj )A j+1 (u) - ££ =1 l k A k (u)} . (49) 

A more convenient expression for such media is given by 

M z {u) = exp[-j z A u (u) du] , (50) 

which applies for piecewise-constant functions A u {uj). It is plausible that 
such an expression applies as well for slowly varying media functions A u {uj). 

As was previously the case, it is useful to separate A u {u) into real and 
imaginary parts, i.e., 

M z (u) = exp[- J z B u (uj) du - i J z q C u (uj) du] . (51) 

Under suitable conditions we may set 

B u (uj) = \a- 1 (u)uj 2 , C u (u) = -v- 1 ^)^ , (52) 
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in which case we are led to 

M z (oj) = exp[— ijga _1 (-u)w 2 du + ij z Q v~ l (u)uo du] . (53) 

Equations (51) and (53) apply to various situations. 

As a simple example of how these equations may be applied let us con- 
sider a signal that penetrates a passive medium with a finite thickness of £ 
and otherwise propagates in free space. Therefore, in this case B u {uo) > 
for < u < £, while for u > £, B u (u) = 0. This example pertains to a 
radar signal that penetrates a canopy of foliage or some other absorbing and 
dispersive medium. Even a comparatively thin, passive medium may be well 
approximated by (53) when A u (uj) is sufficiently large for w^O. Under the 
stated conditions it follows, when z > £, that 

M g (u) = exp{-£cu 2 /2d + i[(z- £)u/c + £u/v}} , (54) 

where c is the speed of light in vacuum, and 

a = £ [j l Q a{u)~ l du]' 1 , 

v = £[j[v(u)- 1 du}' 1 . (55) 

Since v(u) < c, it follows that v < c. For simplicity in our example, let us 
assume that v ~ c. Finally, when £/a is sufficiently large, we are led to 

hit) = VdJ~£ e-^-^l™ J f (s) ds , (56) 

defined for z > £. If F (0) = 0, then appeal can be made to (42). 

This equation should be compared with Eq. (36) that corresponds to 
propagation within a single, homogeneous passive medium. The temporal 
signal in both cases is Gaussian and apart from the overall amplitude, the 
shape is principally determined by the passive medium. This result holds 
whenever the medium acts as a low pass filter for an initial signal that, 
comparatively speaking, has a broad spectral content including u = 0. 

4. Stochastic media 

Under certain circumstances, it is appropriate to regard the passive medium 
as a stochastic variable and to define the output signal as the ensemble av- 
erage of a stochastic output. A description of this kind may arise if the 
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output signal is the superposition of multiple signals from a rapidly chang- 
ing collection of media profiles. More specifically, such a situation may arise 
in the case of a coherent superposition of multiple radar pulse returns that 
have traversed a rapidly fluctuating medium. While details of the fluctuating 
medium may be too difficult to describe, certain overall features - stochas- 
tic invariants - may provide gross characteristics with which selected signal 
properties may be analyzed. 

Let us focus on propagation through a stochastic but otherwise homoge- 
neous medium. In mathematical terms this means that M z (u), and thereby 
A(u), are regarded as suitable stochastic variables. The observed signal f°(t) 
is determined by the expression 

where ((•)) denotes an ensemble average over medium characteristics. Let 
us illustrate the kind of behavior that may occur. For that purpose we again 
suppose that 

A(u) = \a- l u 2 - iv~ x u (58) 

and we assume that a is the only random variable. As an example we suppose 
that the variable a, which can only be positive, is distributed according to a 
probability density given by 

pfa' 1 ) = J^(b a - 1 ) m e- ba ~ 1 , 0<a- 1 <oo, (59) 
ml 

where b is a parameter which essentially sets the scale. In particular, the 
mean value of a -1 is given by 

{a~ 1 ) = r a-Va- 1 ) da" 1 = IH + I . (60) 
Jo o 

As a consequence, we find that 

f o (t) = J_ f -iu,(t-z/v) f oM duj ( 61 \ 

Even before evaluating the Fourier transform it is clear that our stochastic 
example has led to a dramatically enhanced pass band at low frequencies. 
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This result is due, of course, to the nonvanishing probability that arbitrarily 
small values of a^ 1 are present in the ensemble, despite the fact that the 
probability that a^ 1 = is strictly zero for any m > 1. 

As outlined in (42), the leading behavior of (61) is given by 

m=<(t)JMs)ds, (62) 

where 

1 f e~ iult 

m°(t) — — 7- ttt^ — rr du . (63) 

zK ' 2n J (1 + ZLU 2 /b) m+1 v ; 

This expression may be evaluated in the following form: 

1 1 m i 

m°(t) = ^^VW^{J2 C\ m) (VW~z \t - z/v\) } e-VW~*\t-*M , (64) 

1=0 

where the coefficients, for < / < m and < m, are determined by the 
relations 

(~i(m) _ 1 

= (2m -1)!! , 

C[ m) = (2m - 1 - I) Cl m - 1] + C^~ 1] , 1 < I < m - 1 . (65) 



5. Conclusions 

Traditional radar applications involve signal generation, transmission, and 
reception, with little or no spectral modification save for a possible spec- 
tral shift (or dilation) due to relative motion. By contrast, radar signals, 
or electromagnetic waves more generally, that traverse passive media with 
normal absorptive and dispersive properties, typically suffer significant spec- 
tral distortion. For example, over a large spectral domain, an exponential 
attenuation effectively removes any signal energy. An exception to expo- 
nential attenuation occurs at zero frequency, when attenuation is frequently 
absent. On the basis of generally applicable arguments, and under a wide 
set of circumstances, it has been clearly demonstrated in this paper that the 
amplitude of the transmitted signal decays as the inverse square root of the 
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distance traveled rather than exponentially. Such behavior is exactly that 
characteristic of so-called Brillouin precursors [1, 3]. While such phenomena 
are frequently discussed in terms of more advanced mathematical techniques, 
it is demonstrated in the present work that they can be largely understood 
in terms of fairly simple pass band arguments. 1 

The simplicity of the present discussion has enabled us to determine the 
form of signals that have traversed inhomogeneous or stochastic media as 
well. The robustness of the characteristic precursor decay rate has also been 
established for these more general media. 

6. Acknowledgments 

Thanks are expressed to T. Olson and S.V. Shabanov for their interest and 
their comments. Additionally, W. Bomstad is acknowledged for his detailed 
comments on the paper. This work has been partially supported by US 
AFOSR Grant F49620-01-1-0473. 

References 

[1] JACKSON, J.D., 'Classical Electrodynamics' (Wiley & Sons, New York, 
1975) 

[2] HOFFMAN, K., 'Banach Spaces of Analytic Functions' (Prentice Hall, 
London, 1962) 

[3] See, e.g., OUGHSTUN, K.E. and SHERMAN, G.C., 'Propagation of 
Electromagnetic Pulses in a Linear Dispersive Medium with Absorption' 
( Springer- Verlag, Heidelberg, 1994) 

[4] KLAUDER, J.R., PRICE, A.C., DARLINGTON, S., and ALBER- 
SHEIM, W.J., 'The Theory and Design of Chirp Radars', Bell System 
Technical Journal, 1960, 39, (7), pp. 745-808 



1 It is noteworthy that my colleague Tim Olson has devised quite another explanation 
for precursor behavior. See his forthcoming paper "Derivation of Precursors via Finite 
Toeplitz Forms". 



15 



